Vortex quantum dynamics of two dimensional lattice bosons 
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Vortices, which are introduced into a boson su- 
perfluid by rotation or a magnetic field, tend 
to localize in a lattice configuration which coex- 
ists with superfluidity pQ [2J [3]. In two dimen- 
sions a vortex lattice can melt by quantum fluc- 
tuations resulting in a non-superfluid Quantum 
Vortex Liquid (QVL). Present microscopic un- 
derstanding of vortex dynamics of lattice bosons 
is insufficient to predict the actual melting den- 
sity. A missing energy scale, which is difficult 
to obtain perturbatively or semiclassically, is the 
"bare" vortex hopping rate t v on the dual lat- 
tice. Another puzzle is the temperature depen- 
dent Hall conductivity cth(T), which reflects the 
effective vortex Magnus dynamics in the QVL 
phase. In this paper we compute t v and cr H (T) 
by exact diagonalization of finite clusters near 
half filling. Mapping our effective Hamiltonian 
to the Boson Coloumb Liquid simulated by Ref. 
[4], we expect a QVL above a melting density of 
6.5 x 10 -3 vortices per lattice site. The Hall con- 
ductivity near half filling reverses sign in a sharp 
transition accompanied by a vanishing tempera- 
ture scale. At half filling, we show that vortices 
carry spin half degrees of freedom ( 4 v-spins'), as a 
consequence of local non commuting SU(2) sym- 
metries. Our findings could be realized in cold 
atoms, Josephson junction arrays and cuprate su- 
perconductors. 

The model - We consider N b hard core bosons (HCB) 
hopping on a square lattice of unit lattice constant and 
size N = L x L y . An external vector potential A modu- 
lates the hopping amplitude (Josephson energy) t. The 
system is placed on a torus with periodic boundary con- 
ditions, as shown in Fig. [T] 

In the spin- 1 representation of HCB, the angular mo- 
mentum raising and lowering operators S^r create and 
annihilate bosons; the occupation number is n r = S* + \. 
The Hamiltonian we study is a gauged XXZ model, 
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Here rj = ±x, ±y, is the link direction on which the 
lattice gauge field A v is defined. Here we only consider 
the superfluid regime of weak nearest neighbor repulsion 
< V <C t. 

In the absence of external magnetic field, the classical 




FIG. 1: The Gauged Torus. Geometry of HCB Hamilto- 
nian Eq. |l]) which serves to extract its vortex mass and 
Hall conductivity. The torus surface is penetrated by a uni- 
form magnetic field of one flux quantum, and threaded by two 
Aharonov Bohm fluxes = (G x , S y ). For one flux quantum 
there is no translational symmetry on the torus. Red circles 
denote cycles of zero flux, and the vortex center R v (0) is 
localized on the antipodal point to their intersection, the null 
point. 



ground state of 7Y, is a ferromagnet in the XY plane 
with a uniform z-magnetization density m z = n b — 1/2. 
The mean field superfluid stiffness is given by p™ f = 
tn b (l — rib). Consequently, the superfluid transition tem- 
perature, which is proportional to p s , is maximal at half 
filling [5]. 

An important distinction between lattice hard core 
bosons and continuum models is the existence of a charge 
conjugation operator C = exp (in S% ) on the lat- 
tice. C transforms boson "particles" into "holes" rii — > 
(1 — rii), and the Hamiltonian into 

C f W[A, n b ]C = W[-A, 1 - n h ], (2) 

where n b = N b /N is the filling fraction. A consequence 
of ([2| is that the Hall conductivity is antisymmetric in 
n b - 1/2: 

a H (n b ,T) = -a H {l - n 6 , T). (3) 

In terms of vortex motion, this relation implies that be- 
low and above half filling vortices drift in opposite di- 
rections relative to the particle current. Sign reversal 
of Hall conductivity is familiar from tight binding elec- 
trons at half filling on bipartite lattices. Here, however, 
the mechanism of sign reversal is different: it arises from 
the hard-core interactions of bosons and occurs for any 
lattice structure. 
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Vortex effective Hamiltonian - Vortices in a two- 
dimensional superfluid act as charges in (2 + 1)- 
dimensional electrodynamics, where the phonons of the 
superfluid become the photons of the electromagnetic 
theory. The duality between vortices and charges is dis- 
cussed in refs. El US [EE] • The role of speed of 
light is played by the phonon speed of sound, which in our 
model is c = y/2ta/h, where a is the lattice constant. The 
vortex centers R^ are modeled as point charges hopping 
on the dual lattice of plaquette centers, coupled mini- 
mally to the gauge field A^ — + a M , where a M is due 
to the average boson density, generating a dual magnetic 
flux per plaquette of 27rn5, and where A M is a dynamical 
field describing the fluctuations in the boson density and 
current. At low energies, vortex dynamics is described 
by a Harper Hamiltonian plus confining potential, 
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Un(IL) is an effective potential which binds the vortex 
to the 'center of mass' location H v (O x ,Q y ). given by 

nana, 



(5) 



is the number of flux quanta, n a are integers, L a 
are the dimensions of the torus. The phase angles Q a 
are proportional to the solenoidal fluxes of the two el- 
ementary toroidal cycles, and a set of 'null points' at 
the intersections of cycles C XyV is defined by requiring 
§ A • dl = 27rn a /N ( j ) , as shown for = 1 in Fig. 1 
c a 

The vortex CM must lie antipodal to one of the iVj null 
points; we study the simple case = 1. The con- 
stant K is calculated variationally from Eq. ([I]) using 
spin coherent states. Minimizing the energy with respect 
to the position of the vortex centered at R, determines 
the effective potential Un(R)- For V = at we find 
K 39.2 £77,5(1 — rib). We stress that inclusion of the 
potential ?7jv(R) is essential in order to extract the vor- 
tex hopping parameters from our small scale numerical 
calculations. 

For a quantitative quantum theory of vortices we need to 
evaluate the effective hopping t v . Since vortex tunneling 
between lattice sites depend on short range many-body 
correlations, we extract t v from exact numerical diago- 
nalizations of TL on 16 — 20 sites clusters, in the presence 
of a single flux quantum. By tuning t v , we fit the low- 
est three eigenenergies E n and eigenstates \^ n ) of TL to 
the lowest states of the effective Harper Hamiltonian Q. 
This assume that |\£ n ) states correspond to zero point 
fluctuations of the vortex positions, since phonons are 
frozen out by the finite lattice gap of order 2irt/\fN. 

Our results for t v (rib, V/t), for N = 20 fit the analytical 
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FIG. 2: (a) The vorticity (V x j) for the first three doublets 
of the HCB model, Eq.Q, with = 1 and = on a 
4x4 lattice. The uniform background vorticity has been sub- 
tracted, (b) Single particle probability density of the lowest 
three excitations of if v , Eq. (|4| with t v = t. 



approximations, 



t Y (n b ,0) 
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The system parameters were varied in the range \rib- \ \ < 
0.2, and V/t < 0.5. We find that at half filling, t Y varies 
very little between the N — 16 and N — 20 lattices. To 
further test our assignment of t v , we compare the vortic- 
ity density (V x j) of eigenstates of H, to the probability 
density of eigenstates of H v . As shown in Fig. [2j using 
the fitted value of t v we obtain similar distributions for 
both sets of wavefunctions. 

Vortex tunneling - In Q we find that near half filling, 
vortices are as light as bosons, t v ~ t. This implies that 
the vortex tunneling rate between two localized pinning 
potentials of strength V, which are separated by distance 
decays exp onen tially as T ~ Ve~ d l x . The localization 
length A oc y/t/V diverges at weak pinning. This result is 
to be contrasted with weakly interacting continuum Bose 
gas. There, the vortex tunnelling rate between pinning 
sites is much smaller, and decays as a Gaussian T ~ 
e -f n b d ^5 From this comparison, we conclude that at 
half filling = 1/2, the lattice and interactions enhance 
vortex mobility considerably. 

Quantum Melting Transition - Having calculated t v , we 
can write down the effective multi- vortex hamiltonian in 
the thermodynamic limit. We drop Un at large N. By 
Q, at half filling there is dual magnetic flux tt per pla- 
quette. In the magnetic Brillouin zone, there is a two- fold 
degenerate dispersion E^ S1 s =|,|. We later return to 
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explain the origin of this 'v-spin' degeneracy. The vor- 
tex effective mass is M" 1 = d 2 E^/dk 2 = t v a 2 /h 2 . In- 
tegrating out the phonon fluctuations produces an 
instantaneous logarithmic (2D Coulomb) interaction be- 
tween vortices, plus retarded (frequency-dependent) in- 
teractions [TT] . These can be represented by a self energy 
TL ret {uo). Since we are interested in the short wavelength 
fluctuations which are responsible for the quantum melt- 
ing of the vortex lattice, we ignore these retardation ef- 
fects. 

Thus, for half filled bosons and a vortex density n v we 
arrive multivortex Hamiltonian 

-^EN 2 + wret H- ( 7 ) 

i 

The single spin version of 7Y mv , after setting H ret — > 0, 
is the Boson Coloumb Liquid studied by Magro and 
Ceperly (MC) [4 by diffusion Monte-Carlo simulation. 
Their dimensionless parameter which governs the phase 
diagram is r~ 2 = 7rn v aQ. We set their ao = ) 1 ^ 2 
as the microscopic length which matches between their 
model and 7Y mv . MC found that below r s w 12 the boson 
lattice undergoes quantum melting. Using our values of 
t v in Eq. (J6|, the critical r s = 12 translates into a vortex 
melting density of 

n c v r < ^6.5 - 7.9y^ x 10" 3 vortices per site. (8) 

This is a suprisingly low vortex density, which implies 
that a QVL can be created at manageable rotation fre- 
quencies for cold atoms, and moderate magnetic fields for 
Josephson junction arrays and cuprate superconductors. 

Hall Conductance - The temperature-dependent Hall 
conductance of the finite cluster is given by the thermally 
averaged Chern numbers [16] : 



cr H (n 6 ,T) = 




E n (®), and |^ n (0)) are the exact spectrum and eigen- 
states of 0. The results are matched at high tem- 
peratures with the ones obtained by the Kubo formula 
[T7] . A typical Hall conductance as a function of filling 
for = 1 is plotted in Fig. [3] At zero temperature, 
cfh = below half filling., reminiscent of the behavior 
in the continuum ctr oc Nb/N v , which holds irrespective 
of temperatures. However, for HCB <th(T, n^) decreases 
with temperature. Moreover, <th reverses sign at half 
filling, as expected by 

Our results show a striking general feature. We find 
that <jh undergoes a sharp transition between <th > 
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FIG. 3: Hall conductance as a function of boson number 
Nb for Hard Core Bosons, Eq. |l]) on the torus. Temper- 
atures vary in intervals of AT = 0.05£. The jump of the 
zero temperature conductance at half filling, is smoothened 
at finite temperatures. Inset: Hall temperature scale as a 
function of density deviation from half filling. Th is defined 
by a H (T H ) = |cr H (0). 



(on < 0) just below (above) half filling. As the tem- 
perature is lowered, the sign reversal of the Hall con- 
ductance happens across a narrower region around half 
filling. This suggests a singularity in the thermodynamic 
system with a vanishing energy scale. We define Tn(n^) 
by <th(Th) = ^cth(O)- In the inset of Fig.|3j we show that 
Th seems to vanish with \rib — ||, although we cannot yet 
investigate this behavior further in larger systems. 

Spin-^ vortices - Half filling is a special density for TL. 
First, the Hall coefficient vanishes by ([3|, which implies 
that the vortices see no static Magnus field. Second, the 
external magnetic field creates a multitude of doublet de- 
generacies. To be precise, for any odd number N^ of flux 
quanta, there are TV (the system size) distinct values of 
AB fluxes Si where all eigenstates are two-fold degener- 
ate. We have found that these degeneracies are associ- 
ated with non-commuting local symmetry operators 

U a = ^U a CP a [K v ], a = x,y. (10) 

C is the charge conjugation (see Eq. (J3|), and U a is a 
pure gauge transformation. P x ^ is a lattice reflection 
about the x(y) axis passing through the vorticity center 
R v . For N discrete AB fluxes, R v can be placed on each 
one of the lattice positions, where [H, H a ] = 0. These 
symmetries follow from the fact that CP a preserve the 
magnetic field and the AB fluxes. If 0^ are tuned by ^ 
to position R v on a lattice site, H a sends TL to itself upto 
a pure gauge transformation (U a )^ . 

A straightforward, though cumbersome, calculation [17] 



4 



yields the commutation rule 

n y n x = (-i) N *n x w . 



(ii) 



We define the vector II = (IP, H y , II 2 ), where U z = 
2iU x H y . For any odd number of vortices it is easy to 
show using (11) that each of the energy eigenstates is at 



least two-fold degenerate. 

Since II 2 = 3/4, they obey the algebra of spin half oper- 
ators. Thus the doublets reflect the Kramers degeneracy 
expected for an odd number of interacting spin half de- 
grees of freedom which we label v-spins. As shown by 
the form of IT*, the v-spins are attached to the vortex 
positions. The z-direction polarization corresponds to a 
boson charge density wave (CDW) modulation. Varia- 
tional calculation shows the CDW to be exponentially 
localized in the vortex core [18]. Thus v-spin interac- 
tions between different vortices decay exponentially, and 
are very weak in the vortex lattice regime. We cannot 
determine their ordering configuration. If they interact 
ferromagnetically, a CDW order parameter can form. In 
any case, whatever the ordering tendencies, the low ex- 
change energy ensures that v-spins excitations play an 
important role in the low temperature thermodynamics 
and transport coefficients of the multi vortex system in 
both lattice and QVL phases. 

Nature of the QVL - Theoretical treatments of lattice 
bosons have found a myriad of vortex-antivortex con- 
densate (VC) phases at all rational boson filling fractions, 
n b = p/q, due to g-fold degeneracies of the Harper hamil- 
tonian on an infinite lattice [191 HOI HH [22l [23]. VC's 
are in effect, insulating phases where the dual Anderson- 
Higgs mechanism produces a Mott gap [7 . Some of these 
phases result in g-periodic CDWs. Therefore a possible 
candidate for the VC at half filling is the bipartite CDW 
i.e. the antiferromagnetic Ising state of ([I]) at V > t. 
However, the QVL we study, which contains a net den- 
sity of vortices, differs from the proposed VC phases in 
two important respects. 



(i) MC U have found that the liquid phase of H mv Q 
has vanishing condensate fraction. If their results (ignor- 
ing retardation effects) is relevant to the QVL, it should 
differ from a charge-gapped insulator. Whether it is a 
metal is an open possibility. Away from half filling, our 
results for a xy show that the vortices are subject to a 
strong magnetic field, which further suppresses their con- 
densation. At low boson fillings and large vortex density, 

i 'n<f> < 1? there is evidence for fractional quantum hall 
phases [241125] , 

(ii) Away from the commensurate filling = p/q, the 
Hall conductivity is expected to cross zero, and be pro- 
portional to the excess density from p/q. We found that 
the Hall conductance has a very different behavior: it 
has only one abrupt jump between a positive value be- 
low, and negative above half filling. We have computed 
the single flux Hall conductance on a finite lattice, but 
we expect it to represent the macroscopic Hall conductiv- 
ity in the QVL phase where superfluid order parameter 
correlations are short ranged. 

Summary - Vortices of hard core bosons near half fill- 
ing are highly mobile logarithmically interacting charges. 
The vortex lattice is expected to melt into a QVL at 
around 7 x 10 -3 vortices per site. At boson density of 
half filling, the vortices are not subjected to an effective 
Magnus field, but carry local v-spin degrees of freedom, 
which effect the low energy correlations. Away from half 
filling, the Hall conductivity exhibits rapid variation ac- 
companied by a vanishing energy scale. Although the 
issue is far from settled, we present arguments that the 
QVL is not a vortex condensate. 
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